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Based on the calculation of the quasiequilibrium statistical sum 
by means of the functional integration method, we obtained a 
nonequilibrium statistical operator for the electron subsystem of a 
semibounded metal in the framework of the generalized "jellium" 
model in the Gaussian and higher approximations with respect to 
the dynamic electron correlations. This approach allows one to go 
beyond the linear approximation with respect to the gradient of 
the electrochemical potential corresponding to weakly nonequilib- 
rium processes and to obtain generalized transport equations that 
describe nonlinear processes. 



1. Introduction 

Equilibrium characteristics and nonequilibrium pro- 
cesses of diffusion, adsorption, and desorption for spa- 
tially inhomogeneous electron-atom systems are de- 
scribed with the help of various theoretical approaches, 
both available and being developed ones. In particu- 
lar, one widely uses the time-dependent density func- 
tional theory (TDDFT) (lrflo|. During the years of its 
development, the TDDFT has demonstrated significant 
achievements and still extends its limits of application 
though with certain problems 
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14). The basis of 



the TDDFT is the Kolm-Sham density functional theory 
15 4l9| . Another theoretical approach is related to the 
hydrodynamic model of surface plasmons for a spatially 
inhomogeneous electron gas proposed in [20 22] with the 
use of the response theory [23( based on the Boltzmann 



kinetic equation. The quantum statistical theory for 
the description of nonequilibrium processes in "metal- 
adsorbate-gas" systems was developed in works 24-26j| 
using Zubarev's method of nonequilibrium statistical op- 
erator (NSO) [27, 28 1 . In particular, a self-consistent de- 
scription of nonequilibrium processes in the atomic and 
electron subsystems was presented in [24] on the kinetic 
level of the description of electron processes. In the pro- 
cesses of adsorption, desorption, and surface diffusion, a 
metal surface undergoes a reconstruction accompanied 
by a variation of nonequilibrium properties of both elec- 
tron and ion subsystems. In this case, the electrodiffu- 
sion, viscothermal, and electromagnetic properties of the 
electron subsystem change in the field of metal surface 
ions. To study the ion and electron structures of a semi- 
bounded metal, a generalized approach that takes the ef- 
fect of discreteness of the ion subsystem into account and 
is based upon the model of semibounded "jellium" 29. 3C 
was proposed in 
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. 26] . It is worth noting that the in- 
fluence of the discreteness of an ion density on the char- 
acteristics of a semibounded "jellium" was considered 
in fl9l. [31 - 3 3 1 by means of constructing a perturbation 
theory with respect to the electron-ion interaction pseu- 
dopotential. However, the linear response of the elec- 
tron subsystem to the lattice potential did not take the 
effects of inhomogeneity of the electron subs yste m into 
account. The approach described in [25,29|, 3(| allows 
one to model the formation of a surface potential and 
to calculate a large statistical sum for the generalized 
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model in terms of the cumulant averages of the "jellium" 
model. In [25], the generalized "jellium" model served as 
a basis for the statistical description of electrodiffusion 
processes for the electron subsystem of a semibounded 
metal with the use of the NSO method, where the only 
parameter of the reduced description was the nonequi- 
librium average value of the electron density. For such 
a system, the quasiequilibrium statistical sum was cal- 
culated by means of the functional integration method 
for the case of the local electron-ion interaction pseu- 
dopotential of the metal surface. In principle, it allowed 
one to obtain expressions for the nonequilibrium statis- 
tical operator in the Gaussian and higher approxima- 
tions with respect to the dynamic electron correlations. 
In [25( 1 . the nonequilibrium statistical operator and the 
generalized transport equation of inhomogeneous diffu- 
sion were obtained for weakly nonequilibrium processes 
(linear approximation with respect to the gradient of the 
electrochemical potential). The same approximation was 
used to deduce an equation for the "density-density" time 
correlation function that determines the dynamic struc- 
tural factor of the electron subsystem of a semibounded 
metal and to demonstrate the connection of this elec- 
trodiffusion model in the linear approximation with the 
TDDFT 

The given study represents the continuation of work 
(25)]. We will obtain expressions for the nonequilibrium 
statistical operator in the Gaussian and higher approx- 
imations with respect to the dynamic electron correla- 
tions, by calculating the quasiequilibrium statistical sum 
by means of the functional integration method. This ap- 
proach allowed us to go beyond the linear approximation 
with respect to the electrochemical potential. For the 
nonequilibrium statistical operator in the correspond- 
ing approximations, we will obtain generalized transport 
equations for the nonequilibrium average value of the 
electron density for strongly nonequilibrium processes 
for the electron subsystem of a semibounded metal. 
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where the first two terms represent the electron kinetic 
energy and the potential energy of electron-electron in- 
teraction, respectively, the third term stands for the po- 
tential energy of ion-ion interaction, and the last one is 
the energy of electron-ion interaction. Electrons in the 
ion field have the charge e, the mass to, and the coordi- 
nates rj, i = 1, . . . , N. By N- lon , we denote the number 
of metal ions with the charge Ze and the coordinates 
Rj (— 00 < Xj, Yj < +00, Zj zqj %a = const, z = Zq 
is the division plane), j = 1, . . . , JVj on . We suppose that 
ions are immobile in the system volume V = SL, where 
S is the surface area of the semibounded metal, and L 
determines the region of variation of the electron coor- 
dinate normal to the metal surface: z £ (— L/2, +L/2), 
S — y 00, L — y 00. We consider that the system is elec- 
troneutral, i.e. 
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N. 



(2) 



In [25], Hamiltonian |1| was presented in terms of the 
collective variables of the electron subsystem of a semi- 
bounded metal specifying the Hamiltonian of the "jel- 
lium" model reference system: 



H = J2E a (p)ai(p)a a (p) 
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2. Generalized "Jellium" Model. 

Nonequilibrium Statistical Operator 

2.1. Hamiltonian of the system 

Consider an electron-ion system that describes a semi- 
bounded metal with regard for the influence of the dis- 
creteness of the ion subsystem. We present the Hamil- 
tonian of the system in the form 



h 2 N 1 
H = -— VA. + - 
2m ^ 2 
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where the primed sums mean the absence of terms with 
q = due to the electroneutrahty condition ([2]) , v k (q) = 
4ire 2 /(q 2 + k 2 ) and /fe(q) are the three-dimensional 
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Fourier transforms of the Coulomb potential and the lo- 
cal part of pseudopotential ^ : 



7-y"^(q)e' q(r||i_r|l:, ' )+i ' £(z,_2j) , 



\Ti - r., 



SL 



q.A: 



Hamilton operator ((3]). With regard for the chosen ge- 
ometry of the model, the value of (f>(r))* will correspond 
to the mixed Fourier representation (/^(q))'. To find 
p(t) (the solution of the Liouville equation), we employ 
Zubarev's NSO method 27,221 and obtain, in the gen- 
eral case, 



w{ti - Ry) 



Ze 



f(Ti - Ri), 



(4) p{t) 



a e(t'-t) AL N (t'-t) 



Pq(t')dt', 
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u = [X^Y-j), v(q\z) = 27re 2 e-«l z l/(7 is the two- 



dimensional Fourier transform of the Coulomb potential: 

y h 2 p 2 

— = - Y, ^(q|z)e" qr ii , ^(p) = t, he a is the energy 

r b q 2m 

of an electron in the state (p, a), 

N ion 



(5) 



stands for the structural factor of the ion subsystem and 
the Fourier transform of the electron density: 

p k (q)= (ai|e to |a 2 )at i (p)a a3 (p-q) ) (6) 

p,ai,a2 



where (ai|...|a 2 ) = / dz ip* a 1 (z) . . . ip a2 (z). cp a (z) 

and s a are the eigenfunctions and eigenvalues of the 
Schodinger equation, respectively, 

(f a {z) = e a tp a (z), 



2m dz 2 



V(r) = V(z) is the surface potential depending only on 
the electron coordinate normal to the division plane. 

The same way as in [25j], the electrodiffusion processes 
in the formulated model are described choosing the av- 
erage value of the electron density operator as the main 
parameter of the reduced description of nonequilibrium 
processes in the electron subsystem of a semibounded 
metal. It is connected with the corresponding inhomo- 
geneous electric field, 



V-E(r;t) = e( e (r)) t 



(7) 



where ((...))' = Sp(...)p(t), p(t) denotes the nonequi- 
librium statistical operator of the generalized "jellium" 
model that satisfies the Liouville equation with the 



where e — ^ +0 after passing to the thermodynamic 
limit, and iLj^ is the Liouville operator corresponding 
to Hamilton operator ©■ The quantity p q (t) denotes 
the quasiequilibrium statistical operator determined by 
the Gibbs method at fixed values of the parameter of the 
reduced description (/9ft(q)) , and the normalization con- 
dition Svpq (t) = 1 hold. In our case, it has the following 
form 



p q {t) = cxp 



-#(*)- 



^EE^*( q; *)M<l)] 
k q / 



(9) 



where <I>(t) = \nZ(t) is the Massieu-Planck functional, 
and Z(t) is the statistical sum of the quasiequilibrium 
statistical operator, 



Z(t) = Spexp 



\ k q / . 



(10) 

//fe(q; t) = /Ufc(q; t) + e^(q; t) denotes the Fourier trans- 
form of the electron electrochemical potential, /Ltfc(q; t) is 
the Fourier transform of the electron chemical potential, 
and ifk (q; t) is the Fourier transform of the local electric 
potential. The quantity jlk (q; t) is determined from the 
self-consistency condition 

{PM)Y = <P*(q)>* 9 (11) 
and the thermodynamic relations 

S-izPk(q;t) 



(P fc (q)) ( 



(12) 



6S(t) 
SS(t) 



-^Mfe(q; t), 



(13) 



SL 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 2 



181 



P.P. KOSTROBIJ, B.M. MARKOVYCH, A.I. VASYLBNKO et al. 



where S(t) is the Gibbs nonequilibrium entropy, 
S(t) = -Sp(ln p,(t))p,(t) = 



2SL 



fc,q 



& q fe >1 ) 



(17) 



lnZ(t) + p((H) t - ^^^ Mfe (q;t)(p fc (q)}*- 



fe q 



(14) 



fe q 



(epfc(q))* = e(pfe(q))' is the average electron charge den- 
sity. As follows from the structure of the nonequilibrium 
entropy, the transport processes in the system within 
the used model are caused by the gradients of the local 
chemical and electrochemical potentials. 

With regard for the structure of p q (t), the nonequilib- 
rium statistical operator can be presented in the form 

p(t)=e f e <t'-t) e iLN(t'-t) e -S(t') dt , = 



= p q (t)+ / e e(t'-t) e iL N (t'-t) : 





where 



s(t') = in z(t) + ^-^EE A*(q; t)p*(q)) (16) 



where B(q,k;t) = A/i on Sfe(q)a;fc(q) -/ife(q;t), tUfe(q) = 
-Z;/ fc (q) + e/ fe (q), = Ep, Q £; «(p)4(p) a a(p) is thc 
kinetic part of the Hamiltonian of the electron subsys- 
tem. 

Applying the functional integration method and con- 
sidering the "jellium" model as a reference system, Z(t) 
can be written down as 1251 



Z{t) = exp 'Kq|0)| Zjeu AZ(t), 

where 

2jeu = Sp{exp(- J 8fl o )T5i08)} 



(18) 



(19) 



is the statistical sum of the "jellium" model of the elec- 
tron subsystem of a semibounded metal that corresponds 
to the equilibrium state calculated in [29|,|30j; 



Si 08) = exp 



Xp fe (q|/3')p-/e(-q|/?') 



(20) 



is the contribution of the electron interaction, where 
(15) p fc (q|/3') = e/5' H °Pfc(q)e-^°, 



AZ(t) = ^-Sp{exp(-l3H )TS 1 (f3)S 2 (f3;t)} 



= (5 2 (/3;t)) je ii, 



(21) 



is the entropy operator. In order to reveal the structure where (( )) jc n = Spj exp(-/Ji?o)TSi(/3)(...) j; 

of the entropy operator, it is necessary to calculate the JC 
statistical sum Z(t) of the quasiequilibrium statistical S^Cft*) = 
operator. With regard for the structure of Hamiltonian 
([3]), Z{t) can be put down as follows: (25| : 



Z(t) = Sp J exp(-/3(tf - ^^'Kq|0)+ 
I q 



f P 1 

Texp J -J d/3'^^B(q,fc;t)p_ fe (-q;^)[- (22) 
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Using the cumulant representation, AZ(t) can be pre- 
sented in the form 



AZ(t) = exp 



xB(q„,fc n ;t)9Jl_ fel ...._ fe7l (-qi.... - q„) 
where 

9ttfei....fc„(qi--qn) = 

= t n <Tp jkl (q 1 |0),....p fcB (q n |0)>f BU 



(23) 



(24) 



are the cumulant irreducible average values of the elec- 
tron density fluctuations calculated with the help of the 
equilibrium statistical operator of the "jellium" model of 
the electron subsystem of a semibounded metal (29l . l30| . 
In particular, the second cumulant has the structure 

2Kfci,fc 2 (qi>q2) = (pfei(qi)Pte(q2))jeu- 



■ (Pkx (qi)) jell (Pfc 2 (q2)) jell 



(25) 



and is connected with the static structural factor 
5(fei,qi; k 2 ,ci2) = (Pfe 1 (qi)pfc 2 (q2))jell of the electron 
subsystem of a semibounded metal. In the Gaussian 
approximation, we obtain 



AZ G (t) = exp 



qiq 2 fcife 2 



xB(q 2 , fc 2 ;i)9^-fei,-fe 2 (-qi, -q2) 



(26) 



is expressed in terms of the second cumulant of the "jel- 
lium" model of inhomogeneous electron gas (29l . l30| . Ac- 
cording to the definition of s-particle electron distribu- 
tion functions [25,29,30], we obtain the quasiequilibrium 
s-particle electron distribution functions in the form 

F s (r 1 ,...,r n ;t) = F s (r 1 ,...,r„)J' ell x 



xB{q n ,k n ;t)Am ( X, 
where 



-fen 



(-qi,-- 



-q«) 



(27) 
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, (-qi, • • 



-q«J 



= M 



0) 

-fcl,..., — k; 



(-qi,...,-q„)- 



-M- ku ... _fe„(-qi, ■ • ■ , -q«)- 

Relations (|27p link the quasiequilibrium distribution 
functions with the electrochemical potential pk (q; i) 
through the corresponding cumulant averages of the "jel- 
lium" model. In view of the structure of AZ(t) (|26|) . 
In Z(t) can be written as 



2S 



n=l x ' qi ...q„ k\ ...k n 



xB(q n ,fc„;i)97l_ fcl ...._ fcn (-qi.... - q„), 



(28) 



where ln^j e n can be calculated in various approxima- 
tions with respect to the electron correlations (29l. l30l| . 
Based on (f2"5| and (ITBl , we obtain the expression for the 
nonequilibrium statistical operator in the general form: 



^ ^ at 

qi-q™ fei...fe, 



-^^B(q 1 ,k 1 ;t')...B(q n ,k n ;t') 



x9Jt_ fc 
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1 

x j dr^(t')p fc (q)p^ r (t')Mfc(q;^)dt', 



where p k (q) = iL N p k (q) = -kq • Jfc(q), and J fc (q) is 
the Fourier transform of the microscopic electron flux. 
The obtained expression represents a sum of the non- 
dissipative and dissipative parts. The first one corre- 
sponds to the operator p q (t), while the second one is 
described by the terms that contain the time derivatives 
of the functions B(q,k;t') and the microscopic fluxes 
/9fe(q). Moreover, the derivative -^jB{q, k; t') can be pre- 
sented as 

^B(q,k;t') = ^(iV ion 5 fc (qK(q) -Mfe(q;i')) = 



^ (q; ' ) = W* (pt(q)> ' 



dt 

where 

Spkjq; t') 
5( P k(q)Y' 



O fc (q)|p_ fe (-q))" 



(Pfe(q))'j an< i microscopic fluxes p(q) of the electron sub- 
system of a semibounded metal. To make the description 
complete, it is necessary to obtain the transport equa- 

(29) 

tion for (pfc(q))' with the help of the non-equilibrium 
statistical operator ([29]) . Based on the structure of p(i), 
one can state that these equations will be nonlinear. The 
parameters p k (q; t) in these equations should be found 
with the use of the self-consistency conditions (ITU) . The 
obtained expression for the quasiequilibrium statistical 
operator with the Massieu-Planck functional (|28|) allows 
one to find the nonequilibrium statistical operator in the 
corresponding approximations, particularly in the Gaus- 
sian one. 

3. Gaussian Approximation 

Here, we will consider approximation (|26p . in which the 
nonequilibrium statistical operator and the transport 
equation for (pfc(q))* will be obtained. The index G 
in all cases means the description of a function in the 
Gaussian approximation. 

With regard for the structure of (f2l?|) . we obtain the 
entropy operator (|14[) in the form 

(30) M ™,.. ~N 



S^(t)=P—J2'^\0)+lnZ jcll - 



^(q)lp-fe(-q)); 1 = /( (/*(q) - (p fc (q))*') x 



SL 



ki,k 2 qi,q2 



^(Mon5 , fc 1 (qi)wfc 1 (q 2 ) -Mfci (qi;*) 



x (p_ fc (-q;r) - (p_ fc (-q))*') J)*dr, 



(31) 



is the quantum quasiequilibrium correlation function, 
p k (q;T)=p r q (t')p k (q)p- T (t'). 

With regard for (|30p . the time derivative of B(q,k;t') 
can be presented in the form 

| 7 i3(q,fc;0 = (p fc (q)|p_ fc (-q))7 1 ^(q))* = 



(32) 



x(MonSfc 2 (q2)wfc 2 (q 2 ) - ^ 2 (q 2 ;t) ) x 
xaJl_fc I ,_fc 2 (-qi,-q2)+ 



(33) 



In order to eliminate the parameters pk (q; t) from this 
formula, we use the thermodynamic relation (|12|) 



<5$( G )(£) 
S-gzPki^t) 

which yields 



= (M(q)}' 



Thus, the nonequilibrium statistical operator (|29|) is a 
functional of the pair quasiequilibrium correlation func- 
tions ([3T|) . equilibrium correlation functions (|24|) . av- 
erage values of parameters of the reduced description 
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x J2 (Sfc'(q') - M**)) 9Jt_ fc ,,_ fe (- q) -q)5 fe (q), (34) 
where 

5*fe(q) -iVion5 fc (q)w fe (q). 

Defining the function Sfl~J, (— q, ~q') inverse to 
5D?_fc._fc' (— q, — q) by the relation 

]T OT-i„( q ,q'0OT fc ^(^q') - 4, fe '<W 

and using (|34[) . we obtain the Fourier transform of the 
electron electrochemical potential as follows: 

Mfc(q;t)=5 fc (q)- x 



J E E(^( c i0) t ^ 1 (qi)^ 1 ,- fe2 (-qi'- c i2) ; 



fci,fe2 qi.q2 



^ fe2 1 (q 2 )(p fe2 (q 2 )) t +/3U-^^{^.(q)- 



SX 



^( Pfe ,( q ')) t 5 fe - 1 (q')^:^,- fe (-q' ) -q)}p fc (q) 



(37) 



Instead of using the nonequilibrium statistical operator 
in the form (|29[) to obtain the transport equations for 
(Pfc(q)} i we apply the nonequilibrium statistical oper- 
ator with regard for the projection, which allows us to 
eliminate the time derivatives of thermodynamic param- 
eters [H^lii. We obtain 



x E (MqOWMORik-kC-q', -q) 



(35) 



pit) = P&Y 



One can see that the Fourier transform of the elec- 
trochemical potential in the Gaussian approximation is 
expressed in terms of the structural factor of the ion 
subsystem and the Fourier transform of the local part of 
the electron-ion interaction pseudopotential. The time 
dependence is described by the average nonequilibrium 
value of the electron density renormalized through the 
structural factor of the ion subsystem, the pseudopoten- 
tial Wfe(q), and the function 3JtZ fe / _*.(— q', — q) inverse 
to the pair irreducible cumulant average value of the 
electron density fluctuation. Substituting ([33)) into the 
expression for the Massieu-Planck functional, we obtain 

& G \t) = In Z G (t) = ^ E '^1°) + ln Z ^ 



4E E^M^qi; 



x 9JT 



> (-qi,-q 2 )5- 9 1 (q 2 )( Pfc2 (q 2 ))* 



(36) 



In this case, the entropy operator (|33[) takes the form 

^ (G)(<) = /3 5E 'Kq|0) + inz jcll - 



f^TS^l- P q {t'))iL N p{t')dt' , 



(38) 



where 



T,(t,t') = ex P+ <-/(!- P ? (t"))iiivdt' 



denotes the generalized evolution operator with regard 
for projection and P q (t') is the generalized Kawasaki- 
Gunton projection operator, whose structure depends on 
the quasiequilibrium statistical operator p q (t). In our 
case, P q {t) has the form 

m P '= (p«(*)-E^|yT(p fe (q)) t )spp'+ 



EiMrSp(Mq)p') 



fc,q 



5{pM)Y 



(39) 



and the operator properties: P q (t)p(t) = p q (t), 

Pq(t)pq(t) = Pq(t), P,(t)P,(f) = P,(t). In order to cal- 
culate the nonequilibrium statistical operator according 

to ()38j) in the Gaussian approximation for (i), we 
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must find, first of all, the Kawasaki-Gunton projection 
operator. Taking into account that 

pfHt) = exp J - (f}?L £ V(q|0) + bxZ iell - 



fei qi -Q2 



1 



where 



^ G )(fe,q;i) = ^fe(q;*) = &(q) 



(42) 



xan:l i _ fe2 (-q 1 ,-q 2 )5 fc - 1 (q 2 )(p fc2 (q 2 ))* + 



A'.q 



k' ,q ; 



^^^j'^HqlH^^-q, -q) L 



fe'q 



(43) 



the nonequilibrium statistical operator can be written 
down in the form 



and 



*<Pfc(q)>* 



fe',q' 



(40) 



P(*)=P^(*)-E / e £( *'- 4) 7f (*,*') 



x /dr(pf)(i')) r /p(fc,q;t')x 



-(^( q ')) t 5-, 1 (q'))OT:^_ fc (-q' I -q)^7 1 (q)pf ) W, 

we obtain the following expression for the Kawasaki- 
Gunton projection operator: 

pP{t) P ' = { P f\t)+ 



fc,q fc',q 



x2n:^_ fc (-q',-q)5,7 1 (q)(p fe (q)) t 4 G) W)Spp'- 



-EE M<*» - <Mq')W(q')) x 

xa7rJ, ) _ fc (-q',-q)5fc 1 (q3Sp(p*(q)//)^ G) (*). 
With regard for (|4"Tj) and the relation 

iL ivP? (t) = ^^ G )(A;,q;t)x 

A'.q 



x(pf ) (0) 1 ~ T W< G >(*,q;f , )< 



(41) 



'9 

where 



I p {k^,t') = {\-V^\t l ))iL N p k {q) 



(44) 



(45) 



is generalized diffusion flow, V^ G '(t) is projection oper- 
ator acting on the operator 

v^ G \t)A^Y. E ^*'(q';*)x (46) 

A/ ,q' fc" ,q" 

xa7l- fc , ) _ fc /,(-q / ,-q")^(q")<P*"(q")^)G. 

where <5p*(q';t) = p*(q') - (^(q'))* (^'^V). 

((— ))(G) = Sp(.../9g (*')) is the averaging with the 
quasiequilibrium statistical operator in the Gaussian 
approximation. In its structure, the nonequilibrium 
statistical operator is a functional of the microscopic 
fluxes Pfc(q), the observable quantities (pfc(q))*, and the 
quasiequilibrium and equilibrium correlation functions 
of the electron subsystem of a semibounded metal. With 
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its help, we obtain the transport equation for (pfe(q))* xS(qa, k$; t)DJl-k 1 ,-k 2 ,-k s (— qi> — Q2j — ^3)+ 
in the form 

d 

^<Pfc(q)}* = (pfc(q)) = / 1 „ x 

9t ^E^^)/^))- (49) 

* In order to eliminate the parameters flk (q; t) from this 

= - E / eE(t ' ~ t] d j G j ( fc > q; fc ': q'; * 5 t'W {G) (k\ q; t')dt', formula, the thermodynamic relation ^ will be applied 

k'tfjfn once again: 

(47) SHt) . . .. t 

x p - = (p*(<0> • 

where o-^fx^t) 

q' k' q''i i') = From here, we derive the equation for £ijt 4 (qi;i): 



</ p (fc,q)T g G (t,t')/p(fc',q';r)); 



0%(q)>' = ^EE (^(qO-zMqi;*))* 



(G) 



fcq.(j fc (q)T g G (t,t')J fc '(q';r))* • qV 

\ / (G) 





^qi,q2 fei,fe 2 



xS* fe2 (q 2 )97t_ fel _ fe2 (-qi, -q 2 )- 



~m ~^\~qt) E E (^i(qi)-Mi(qi;*) 

= fcq. J D, ( / G / ) (fc,q;fc',q';i,0-q'fc / , (48) 31 V 5i / qi ^ q3 fcl ^ 3 V 

Dj j (fc, q; A; , q ; t, t ) is the generalized diffusion coeffi- 
cient of electrons in a semibounded metal calculated us- x (Sk 2 { c i2) — /2fe 2 (q 2 ;i)J 5fc 3 (q3)x 
ing the quasiequilibrium statistical operator in the Gaus- 
sian approximation. 



4. Approximation -Bfe(q; f)J3fc' (q'; t)Bk" (q"; i) 



x97l_ fcl _ fc2 _ fc3 (-qi, -q 2 ,-q 3 )- (50) 

T , .1,1 • n , 1 ^, In its structure, this equation is square with respect to 

Let us consider the next approximation after the Gaus- .„ . _ , . J_ . . 

r . , . . i ,1 the functions /1 i. q; t). lo solve it approximately, we re- 

sian one lor the quasiequilibrium statistical sum or the \ . , , 

Massien-Planck functional Tn this case we obtain plaCe 0ne ° f the <l uantlties Mfcifai;*) on the right-hand 

IVldtaOlCll L lailCA lUIlCLlOIiai ll^OIJ . Ill LI11S Last. Wt OULaill • i r l l • c l* if l • l /~< 

, . r . . r 77"^ side of the quadratic form by its value found m the Gaus- 
the following expression tor the entropy operator: 

25 



sian approximation (|35p . Then we obtain a linear equa- 



N , tion for the function ui.(a:t): 

«'(*) = feE Kq|o) + in^ jeU - 



(MW—^E^q)-^';*)) 



E E ^(qa»*i;*)S(qa,fe;t): 



qi,q2 fel,/C2 



xwt-kx ,-fc 2 (-qi> -qa) 



fc',q 



xG_ fe '-fe(-q,-q;t), (51) 
where 

G_ fcl ,-fc 2 (-qi, -qa;*) = 



•J- V^^/ qi ,q 2 ,q 3 fel ,fe 2 ,fc 3 = S , fc 2 (q 2 )Wt-fe 1 -fe 2 (-qi,-q 2 )- 
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w E E w(q'))*^ 1 (q')x x (p-wr <p»'i<nr co 1 + 

q'.qs fc'ifca 



xOTI^ _ fe f-q',-q 2 )x 



xS , /c 3 (q3)OT-/ Cl ,-.fc 2 ,_fc3(-qi,-q 2 ,-q3). (52) 



\ fc,q \ v K ,q 



xG«(q,q';t) Pfc(q) , (54) 



As is seen, the function G—k 1 ,—k 2 (~ c lU~ < l2]t) depends / / 

on time through the observable quantities (pi c i(q)) t . It w here 
also depends on the structural factor of the ion subsys- 
tem Sk'(q), the Fourier transform of the local part of G^/ fe „(q',q";t) = S~] G^/^fq qi;t)x 
the electron- ion interaction pseudopotential Wfe(q), and fei,fe 2 qi,q2 
the cumulant irreducible average values of the electron 

density fluctuations: pair 97l_fc',_fc 2 (-q, -q 2 ) and triple , 
ones9Jt_ fel ,_fe 2 ,-fc 3 (-qi,-q2,-q3). The second term on XJ ' I -fci.-*iil- < li>- < »; ( 'te,fc'H<l3i<! 
the right-hand side of Eq. (|52p involves the renormaliza- 

tion of the triple electron correlations through the pair ^(3) / / // '".-A — 
ones that make the dominant contribution in the Gaus- k f M" ,k"' ^ > ^ ' ^ » 
sian approximation (see the previous section). 

Defining GlJ. i _ fe2 (-qi, -q 2 ;i) as a function inverse _ q-i (q' ) q 1 ^w 

to G_ fa ,_ fa (-qi',-cte;t) by the relation . >,. i; „ ,,,,,, ^ 

E G =L - fe "(-qi,-q";*)G- fe ",- fc2 (-q",-q 2 ;t) = 

k",q" ' xG^ 2 (q",q 2 ;t)G^, fa (q'",q 3 ;t)x 

= 4i,fc2 5 qi,q2 xOT-fc! _fe 2 ,_fc s (-qi, -q 2 , -qa)- 

and using Eq.©, we derive the following expression Thege functions involve the dynamic renormalization 

for the Fourier transform of the electron electrochemical of the cumulant i rre ducible average values of the elec- 

potential: tron density fluctuations: pair SJl-fe'^^^ (—q, —q 2 ) and 

jSfc(q;t) = Sk(q)- triple ° neS 2,t-kl -fe,-fe 3 (-qi: -q2, -Qs) through func- 

tions (j52")) . With regard for the entropy operator 
the quasiequilibrium statistical operator reads 

E^^^'Glt^C-q'.-qjt). (53) p £W) (t) = exp ( _ £ V(q | 0) + ln 

fc',q \ q 

Now, with regard for (p)3")) , the entropy operator P^l) can 

be presented in the form 1 ; / w-t ^vf2l , / /; ,\ 

-jL L (Pfe'(q)) G fc''fc"(q ,q 

£'(t)=/3g $>(#)+ In Z jell - fe ' fe " q ' q " 

9 



X^q")}^ 



4 E E^^)*^!*"^''^*)^'^'))^ 



3! 

k' ,k" ,k"' q ,q",q'" 



+ 4 E E G^„, fc ,„(q',q",q'";i)x 

3 ^w»' q w' x(p fe Kq')) t (M"(q // ))*^"'(q / ")) i + 
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; (^-^E(^( q )-(i)"E^( q ')) 

\ fc,q k' ,q ; 



statistical operator in approximation (1551) . Moreover, 



P q (i') has the following structure: 

< +l w = (V +1) w- 



xGfc4(q,q';t))p fc (q))J, (55) 

where the index "(G + 1)" denotes the third order with — ^ J (q', r)G,r,).,(q, q'; t)- 

respect to the observable parameters in the quasiequilib- k,q I 
rium statistical operator. In addition to the pair (Gaus- 
sian) one, it also allows for the cubic dependence on the ^ _ 
parameters of the reduced description (pfc(q))' with dy- — (Pfe'lQ )) G k , k (q_ 
namic renormalizations in the functions G k , k ,, (q', q"; i) 



i(3) / / // „lll. 



fc,q 



and Gi/ k „ k ,„(q' ,q" ,q"';t). Since * / / iwt, , nwt 

X 2L L^'W'^'W" 

k' ,k" q'.q" 

xG^, fc (q',q",q;t)}(p fe (q))Vf +1) W)Sp(p')x 



(/3^(q',r)G fe -, 1 fc (q,q';t)- 



x y(pf +1) ) T (i)p fc (q)(pf +1) ) 1 - T (t)dr, (56) 
o , 

where 

<G + iu ^ /9 , , P \- , , -(^(q')) t G^ fc ( q ',q;t))x 
^ G+1 )(fc,q;i) = ^/2fe(q;t) = W S fc (q)- 



SL^™ ' SL 



■§ij Ew^'^k-kH'-^*)}' (57) 

fe'q J 



we obtain the following expression for the nonequilib- xG k> ,k" A^' ^" i^ 1 ) j S v(Pk(<l)p') P g G+1 (*)• (59) 
rium statistical operator in approximation ([55| : 



As compared to the result of action of the operator 

t { C} 

ft')/ i W m ^ ne Gaussian approximation, it already in- 

eludes the third order with respect to the parameters of 
fc ii-oo the reduced description. With the help of the nonequi- 

librium statistical operator (|55p , we obtain the transport 
equation for (/9fc(q))* in the form 



xW {G+1) (k,q;t')dt', (58) 



t 

fe '.q'-oo 



where the Kawasaki-Gunton projection operator 
Pq G+1 \t') and the respective evolution operator 

T^ G+1 \t,t') are calculated with the quasiequilibrium x W {G+1 \k', q; t')dt'+ 
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t 

+ E / e £ (*'-*)(p fc (q)Tf)(t,OPf +1) (^)x 



X^(q'5 T)Yl G+1) W^ G+1 \k\ q'; df , (60) 
where 

Dfj +1) (k,q; k',q';t,t') = 

= (p k (q)Tf(t,t>)p k ,(q',T)) t ( G+1) = 

= k q ( J fc (q)T g G (t , t') J* (q ';r) ) f G+1) q' fc' = 



= fc qJ D^ +1) (fc,q;fc',q';t,Oq'fc', (61) 

Dj^ +1 ^ (k, q; fc', q'; i, t') is the generalized diffusion coef- 
ficient of electrons in a semibounded metal calculated 
with the quasiequilibrium statistical operator in approx- 
imation (|55p . With regard for the structure of functions 
(|52|) and ([57)1 and the Kawasaki-Gunton projection op- 
erator (|5^| . we can conclude that Eq. (fBT))) is nonlinear 
with respect to (^(q))*. 

5. Conclusions 

Electrodiffusion processes in the electron subsystem of a 
semibounded metal are described on the basis of the gen- 
eralized "jellium" model with the use of the NSO method, 
where the only parameter of the reduced description is 
the nonequilibrium average value of the electron den- 
sity. Applying the functional integration technique, we 
have calculated the quasiequilibrium statistical sum for 
such a system in the case of the local pseudopotential 
of electron-ion interaction in a metal in the Gaussian 
and higher approximations with respect to the dynamic 
electron correlations. They are used to obtain expres- 
sions for the nonequilibrium statistical operator in the 
Gaussian and higher approximations with respect to the 
dynamic electron correlations, which makes it possible to 
go beyond the linear approximation with respect to the 
electrochemical potential. In the respective approxima- 
tions for the nonequilibrium statistical operator, we have 
derived the generalized transport equations (generalized 
diffusion equations) for the nonequilibrium average value 
of the electron density that can be applied to the descrip- 
tion of strongly nonequilibrium processes for the electron 



subsystem of a semibounded metal. The generalized dif- 
fusion coefficients for electrons in a semibounded metal 
that enter the corresponding transport equations are cal- 
culated with the quasiequilibrium statistical operator in 
the respective approximations: Gaussian one (|40|) and 
approximation (|55l) . An important point in such an ap- 
proach is that the time correlation functions and the 
generalized diffusion coefficients are calculated with the 
quasiequilibrium statistical operator in the correspond- 
ing approximation and represent functionals of the ob- 
servable quantities (/?fe(q))' of a certain order. Of special 
interest in this approach are the investigations of the dy- 
namic structural factor for the nonequilibrium electron 
subsystem of a semibounded metal. 
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/JO CTATHCTHTHOrO OIIHCy ETEKTPOflH©y3IHHI!X 
nPOLIECIB EJIEKTPOHHOI m/JCHCTEMH 
HAIIIBOBME>KEHOrO METAJIY 
B y3ArAJ!BHEHIH MO/JEJII ">KEJIE" 

77.77. Kocmpo6iu, B.M. MapKoeun, A.I. BacuAeHKO, 
M.B. ToKapuyK 

P e 3 K) m e 

3a flonoivioroio MeTOfly <pyHKn,ioHajibHoro iHTerpyBaHFra OTpiiMaHO 
HepiBHOBajKHHii CTaracTHHHHH onepaTop flJIH ejieKTpOHHOl niflCH- 
ctgmh HaniBo6MejKeHoro MeTajiy b y3arajibHeHift MOflejii "scejie" y 

raycoBOMy Ta Biinpx Ha6jiH»ceHHax 3a flMHainnnHMH ejieKTpoH- 

roiMii KopejiHiriHMii npn po3paxyHKy KBasipiBHOBajKHoi' CTaTHCTii- 

hhoi cyMH. TaKHii niflxifl flae MOjKjiHBiCTb bhhth 3a MejKi jiiHiii- 

Horo Ha6jiH»ceHH5i 3a rpaflieHTOM ejieKTpoxiiviiHHoro noTeHn,iajiy, 

aKe Bi/i,noBi/i;ae cjia6o HepiBHOBajKHHM npon,ecaM, Ta OTpiiMaTii 

y3arajibHeHi pibhhhhh nepeHOcy, hki onncyiOTb HejiiHiiiHi 

npoi;ecH. 
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